Cold Atom Optical Lattices as Quantum Analog Simulators for Aperiodic 
One-Dimensional Localization Without Disorder 
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Cold atom optical lattices allow for the study of quantum localization and mobility edges in 
a disorder-free environment. We predict the existence of an Anderson- like insulator with sharp 
mobility edges in a one-dimensional nearly-periodic optical lattice. We show that the mobility edge 
manifests itself as the early onset of pinning in center of mass dipole oscillations in the presence of 
a magnetic trap which should be observable in optical lattices. 
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Optical lattices incorporating ultracold atomic conden- 
sates are rapidly becoming ideal quantum systems for 
studying various model Hamiltonians developed earlier 
for studying solid-state phenomena. This is primarily 
due to the extraordinary level of precise tunability that 
experimentalists have achieved in controlling the param- 
eters (e.g. hopping, interaction and lattice periodicity) 
of the optical lattice, which makes it possible for the 
cold atom optical lattice to operate as an ideal quan- 
tum analog simulator for various many-body condensed 
matter Hamiltonians. By contrast, ideal model Hamilto- 
nians (e.g. Hubbard and Anderson models) often poorly 
describe solid-state systems since experimental control 
over complex condensed matter systems is, in general, 
quite limited. In addition solid-state systems are invari- 
ably contaminated by unknown disorder, defects, and 
impurities whose effects are not easy to incorporate in 
model Hamiltonians. The cold atom optical lattices are 
therefore becoming increasingly important in advancing 
our knowledge about the quantum phase diagram and 
crossover in model many-body Hamiltonians of intrinsic 
interest. Examples include: the Bose-Hubbard model [l| , 
the Tonks-Girardeau gas Q , and the BEC-BCS crossover 



In addition to studying strong correlation effects (e.g. 
the superfiuid-Mott insulator transition in the Bosc- 
Hubbard model) in many-body Hamiltonians, cold atom 
optical lattices also offer ideal systems for studying quan- 
tum transport phenomena including ballistic quantum 
transport [100 and quantum localization 0, la, 0, El ■ 
The latter may be more generally classified as metal- 
insulator transition phenomena with a direct relationship 
to the solid-state. The distinction between a "metal" (i.e. 
a system with finite resistivity at zero temperature) and 
an "insulator" (i.e. a system with infinite zero tempera- 
ture resistivity) is purely quantum. Broadly speaking, 
there are four classes of metal-insulator transitions in 
quantum lattice systems: Metal-band insulator transi- 
tion in an ordered periodic lattice arising from the chem- 
ical potential moving into energy band gaps; interaction 
induced metal-insulator transition as in the Mott transi- 



tion; disorder induced quantum localization (i.e. Ander- 
son localization [ljj); and quantum localization in aperi- 
odic (but deterministic) potentials in disorder-free lattice 
systems. 

In this paper, we establish that very general as- 
pects of the metal-insulator transition phenomena (in 
the disorder-free environment) can be directly experi- 
mentally studied in aperiodic cold atom optical lattices 
with the tuning of experimental parameters leading to 
the observation of both band and quantum (Anderson- 
like) localization in the same system but in different pa- 
rameter regimes. Such an experimental study of localiza- 
tion or insulating transitions in deterministic aperiodic 
systems is impossible in solid state lattice systems since 
disorder (which leads to direct Anderson localization) is 
invariably present in solid state systems overwhelming 
any subtle localization effects arising from deterministic 
aperiodic potentials. In particular, all states are local- 
ized in one-dimensional systems in the presence of any 
disorder whereas one-dimensional aperiodic potentials al- 
low for the existence of extended quantum eigenstates. 
This makes one-dimensional optical lattice systems par- 
ticularly interesting from the perspective of localization 
studies in deterministic aperiodic potentials since such 
studies in the corresponding one-dimensional solid-state 
systems are essentially impossible due to disorder effects. 
We therefore consider aperiodic quantum localization in 
one-dimensional optical lattices, conclusively establish- 
ing the feasibility of studying this unusual phenomenon 
in cold atom optical lattices. 

The single-particle quantum localization problem in a 
deterministic quasiperiodic potential (i.e. two lattice po- 
tentials with mutually incommensurate periods) has a 
long history 0,0]. In particular, localization properties 
have been extensively studied in the Harper (or, equiv- 
alently, Aubry) model which has an intriguing self-dual 
point where the eigenstates form a multifractal Cantor 
set spectra and are neither localized nor extended. Away 
from the dual point conventional wisdom dictates that 
all states, as a function of the chemical potential, are 
either all extended or all localized, depending on the mu- 
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tual strengths of the potential and hopping terms. Such 
Harper model type quasiperiodic potentials therefore do 
not allow for the existence of a mobility edge separating 
extended states (above the mobility edge) from localized 
states (below the mobility edge) which is the hallmark of 
the Anderson localization transition in three-dimensional 
disordered system. Central to our work is the conclusive 
theoretical demonstration of a class of one-dimensional 
optical lattice systems where the deterministic lattice po- 
tential does allow for the existence of a mobility edge in 
one dimension [l^j . which cannot happen through An- 
derson localization with disorder. This class of models 
distinguishes itself from other models discussed in the 
context of optical lattices through the formation of a 
metal-insulator mobility edge rather than a metal-band 
edge. We find that: 1) Direct numerical simulation and 
an analytic WKB approximation provide conclusive evi- 
dence for a rare metal-insulator mobility edge in a one- 
dimensional model, the nearly periodic Harper model. 
2) Transport measurements in suitably designed, one- 
dimensional optical lattices can exhibit the mobility edge. 

We consider spinless fermions (or equivalently hardcore 
bosons sufficiently near the Tonks-Girardeau regime) in 
the lowest band of a one-dimensional, tight binding lat- 
tice with external potentials: 



Un+l 



+ {V n + V D f n + nn 2 )u n = Eu n , (1) 



where the amplitudes u n multiply the Wannier states 
at sites n in the real space wavefunction ^(x) — 
J2 n u nw(x — n). We work in units of the hopping matrix 
element, t = 1, and lattice spacing of the primary lattice 
defining the tight binding problem, a = 1, unless other- 
wise noted. The statistics of spinless fermions implicitly 
allow for an arbitrary on-site interaction in the above 
single-band model. In the absence of an external poten- 
tial the solutions form extended states, u n = uq exp(in0), 
with band energies E = 2cos(0), for < <\> < it. The 
band edges lie at E = ±2. In the presence of an os- 
cillatory modulation of strength V, much weaker than 
the primary lattice, we can ignore modifications to the 
hopping. In this limit we impose a secondary lattice: 
V n = ycos(27ran) — V. For a irrational the additional 
potential establishes an incommensurate pseudorandom 
model, the Harper model (for O = and Vd = 0). The 
potential Vofn adds disorder where f n is a random num- 
ber satisfying < /„ < 1 for each site. The confinement 
potential, fin 2 , applies to optical lattice systems. 

According to the Aubry-Andre conjecture 0] the 
Harper model exhibits a metal-insulator transition at the 
self-dual point V = 2. For V < 2 all states are extended 
while for V > 2 all states localize (the states at V = 2 
are critical). The localized states are characterized by a 
nonzero Lyapunov exponent (inverse localization length), 
j(E), where u n (E) ~ exp(— 771), and gaps in the energy 
spectra. While exceptions to the Aubry-Andre conjec- 
ture have been rigorously proven for specific values of a 




FIG. 1: Top panel: Lyapunov exponent versus energy for 
the nearly-periodic Harper model, Eq. 10 with V = 0.5, 
e — 0.005, and fi = 0. The solid (dot-dashed) line shows 
the disorder-free, Vd = (disordered, Vd = 1) case. Bot- 
tom panel: Level number versus energy of the disorder-free, 
nearly-periodic Harper model for two characteristic values of 
e with V = 0.5, N = 3000, and fi = 0. The lower curve is 
shifted downward by 100 levels for clarity. 



|l7| . we discuss here an additional and experimentally 
relevant counter example defined by: a — m ± e, for 
integer m and 



N- 1 < e « 1, 



(2) 



with N sites. In the limit N — ► 00 the secondary lat- 
tice defines a slowly varying, nearly-periodic potential. A 
similar, slowly varying potential has been considered in 
the context of one-dimensional localization in quasiperi- 
odic systems [l8j ]. 

In the limit defined by Eq. J5J the eigenstates of Eq. (JTJ 
with = and Vd — display Anderson-like localiza- 
tion where we expect to find only extended states. To 
see this consider 7 defined in the limit, N — ► 00 fl9| : 
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The first equality allows us to use the transfer matrix 
method to calculate 7 for large system sizes. The solid 
line in the top panel of Fig. ^ plots the Lyapunov ex- 
ponent versus energy for N — 10 7 and V = 0.5. The 
additional potential, V n , shifts the lower band edge to 
E = — 2 — 2V while leaving the upper band edge at 
E = 2. We see extended states in the center of the band, 
— 2 < E < 2 — 2V, with 7 = 0, as expected from the 
Aubry-Andre conjecture. However, near the band edges, 
-2-2F < E < -2 and 2-2V < E < 2, the states local- 
ize, 7 > 0. The points E — —2 and 2 — 2V define mobility 
edges which are unexpected in one dimension but found 
in three-dimensional models with disorder. The localiza- 
tion is, in this sense, Anderson-like. We find that, for 
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N = 10 7 , the mobility edges persist for rational and irra- 
tional values of e from 10 -5 to 10~ 2 . We conjecture that 
in the limit Eq. irrational numbers are approximated 
by rational numbers up to a number much smaller than 
TV -1 . For N — ► oo, the spectra can contain an infinite 
number of infinitely small gaps and therefore localized 
states eliminating the distinction between an incommen- 
surate and commensurate system 0, ^) . 

The unexpected insulating behavior coincides with a 
devil's staircase-like structure in part of the energy spec- 
trum 0, Q . The second equality in Eq. shows that 
a degeneracy at Ej supports non-zero ■y(Ej). The lower 
panel of Fig. ^ plots the level number as a function of 
energy as determined by exact diagonalization of Eq. Q 
for N = 3000, V = 0.5, V D = 0, and fl = 0. The top line 
(e = 0.005) shows smaller gaps and narrower steps than 
the bottom line (e = 0.05) suggesting that the localized 
states develop a gapless insulator in the limits e — > and 
e»JV- x . 

We can understand the insulating states in a "semi- 
classical" approximation where e plays the role of H. We 
analyze the behavior of each regime as a function of en- 
ergy. At low energies, E < —2, the slowly varying po- 
tential confines low energy states near the potential min- 
ima defined by V n . Very little tunneling between min- 
ima forces localization. Intermediate energies,— 2 < E < 
2 — 2V, see a smaller barrier between minima allowing 
for extended states and, therefore, the first mobility edge 
at E = — 2. A second mobility edge forms at E = 2 — 2V 
when states localize at the secondary lattice maxima. At 
first this seems counterintuitive but can be understood 
in a WKB approximation based on the slowly varying 
nature of V n . A similar analysis was performed for a dif- 
ferent model in Ref. ■ Our results show that the high 
energy states, E > 2 — 2V, moving energetically above 
the lattice slow when passing secondary lattice maxima 
to force localization. We have checked that our analysis 
based on the WKB approximation reproduces the solid 
line in the upper panel of Fig. ^ As an additional check 
we can, in a continuum approximation jl2l |20| . define 
a position variable, n — > x, and a difference operator, 
u n +i + Ufi-i ~~ * 2cos(p)w(i) (with p = id/dx), to give 
the semi-classical Hamiltonian: Hcl — — 2cos(p) + V(x), 
with the replacement V n — > V(x). The phase trajectories 
of Hql produce extended and localized states (and there- 
fore mobility edges) in the regimes obtained in Fig. ^ 

We now discuss the possibility of observing this unique 
type of localization. In the solid state a necessary cor- 
rection to the Harper model includes disorder where we 
add to V n a potential of the form: Vb/n- For V — (and 
fl = 0) this defines the one-dimensional Anderson model 
where we expect all states to localize for arbitrary Vn- 
However, for V ^ 0, the states (otherwise extended in 
the Vb = case) have a small localization length which 
could allow some remnant of a mobility edge. The dot- 
dashed line in the upper panel of Fig. ^ plots 7 for the 



same parameters as the solid line but with Vb = 1.0. We 
find that a finite amount of disorder obscures the position 
of the remnant-mobility edges while localizing all states. 

In what follows we consider an essentially disorder-free 
manifestation of Eq. Q : one-dimensional, cold atom op- 
tical lattices. The interference of appropriately detuned 
lasers of wavelength A = 2a can give rise to our tight 
binding lattice with a sufficiently strong lattice height Vz, . 
To create a secondary modulating potential, V n , consider 
an additional pair of lasers at angles 9 and tt— 9 to the pri- 
mary lattice with wavelength A' and amplitude V[. The 
additional lasers interfere to modulate the energy of the 
nth site by: V^cos(27ma) \w(u)\ 2 cos(27rua)du — V[, 
where a — (A/A')cos(0). For small angles we can retrieve, 
up to an overall constant, our nearly-periodic Harper 
model with m = A/A' an integer and e ps —\9 2 /2X'. For 
realistic parameters: Vl = 5-Er, V[ = 0.1Er,9 = 5°, 
and A = A' (where Er is the photon recoil energy), we 
find t fa 0.065£ fl , V « 0.055E R , and e ps 0.004 yield- 
ing the appropriate parameter regime. Furthermore, we 
find that, in the limit of Eq. J5J, fluctuations in the rela- 
tive phase do not alter the position of the mobility edge. 
We now include an important modification to the model 
which accounts for realistic finite size effects. 

A crucial addition to the Harper model in optical lat- 
tices is the parabolic confinement: fin 2 , which leads to 
a finite particle number. We find that weak confine- 
ment leaves the mobility edges intact. To see this con- 
sider the local Lyapunov exponent: j L (Ej) — (2Nql + 
X^=-at cl m \ u n+i/u n \, where the semiclassical lim- 
its of the parabolic trap define the number of states par- 
ticipating in transport, 2Nql + 1- The classical turn- 
ing points give A^cl = 2|xcl(-£')| and Eq. J5J) becomes: 
(2Ncl)^ 1 « e « 1. To determine xcl we set p = in 
H CL with V{x) = Vcos(27rai) - V + fix 2 . For fl - 10~ 5 
we find 2A*cl ~ 10 3 . In the limit fl — > we retrieve the 
usual Lyapunov exponent, j L — > 7. Fig.[3plots the local 
Lyapunov exponent as a function of energy for N = 10 7 , 
e = 0.005, V = 0.5, V D = 0, and fl = 10~ 5 . The mobil- 
ity edges remain even with a reduced number of states 
comprising the system. The inset shows the normalized 
density profile as a function of site number for three dif- 
ferent chemical potentials, (i. At zero temperature we 
include states with E < fx. For fi = 0.5 (dashed-dotted 
line) we find extended states with some modulation due 
to V n . For fj, = 1.5 (dashed line) we have crossed the 
mobility edge and the density pins to unity at some lat- 
tice sites. Here the formation of a mesoscopic version of 
the Anderson-like gapless insulator fixes the density. For 
fx = 3.0 we enter the band insulator regime which fixes a 
large fraction of the states at integer density. 

Dipolc oscillations in harmonically confined atomic 
gases serve as a direct probe of localization ■ A small 
shift in the center of mass results in harmonic oscillations 
in the absence of an external lattice. The presence of one 
or more weak lattices allows for weakly localized states 
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FIG. 2: The local Lyapunov exponent versus energy for the 
same parameters as the upper panel of Fig. Q but with an 
additional parabolic confinement Q — 10 -5 and no disorder, 
Vd = 0. The inset shows the normalized density of the same 
system as a function of lattice number for several different 
chemical potentials. 



which can suppress oscillations and lead to an effective 
under-damping of the center of mass motion. The ad- 
dition of strongly localized states can, in the absence of 
dissipation, eventually pin the center of mass to effec- 
tively over-damp the center of mass oscillations. Strong 
experimental and theoretical evidence supports the pos- 
sibility that band localization has indeed been observed 
in fermionic, one-dimensional optical lattices 0. Similar 
evidence also suggests such behavior for strongly inter- 
acting bosons |8|. 

We now study the onset of the gapless Anderson-like 
insulator and its effect on center of mass oscillations. 
Consider the center of mass to be displaced A lattice 
sites at some initial time T = 0. For extended states, 
the center of mass position, X(T), averages to zero for 
long times while localized states should pin the center of 
mass position, X ~ A. The center of mass position can, 
for some parameters, demonstrate complex, damping-like 
behavior as function of time making a damping constant 
ill-defined. To extract a simple quantity to be compared 
with experiment we calculate the long time average of 
the center of mass position, < X >oo, as a function 
of chemical potential by diagonalizing Eq. with a 
parabolic potential, n = 1(T 5 , for N = 3000, A = -3, 
and Vd — 0. As an intermediate step we require degener- 
ate eigenstates (localized at the edges) to simultaneously 
diagonalizc the parity operator since our system possess 
reflection symmetry about the origin. The dashed line in 
Fig. |21 plots < X >oo as a function of chemical potential 
in the absence of a secondary lattice, V = 0. For fi < 2 
the extended states perform several oscillations about the 
trap center but over long times average to zero displace- 
ment. Above the band edge (labeled B.E.), for /i > 2, 
localized states near the edge pin the center of mass near 
A. For fx > 3 the system never leaves its initial position. 

A second weaker lattice causes a mobility edge to form 
energetically below the band edge. The solid line in Fig. [3] 
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FIG. 3: The long time average of the center of mass position 
as a function of chemical potential after an initial displace- 
ment of three sites, A = —3, with no disorder, Vd = 0. The 
dashed line is calculated from the bare tight binding model 
with no secondary lattice, V = 0, and Q, — 10 -5 . B.E la- 
bels the upper band edge. The solid line is calculated for the 
same parameters but with V = 0.5. The early onset of the 
gapless insulator occurs at the upper mobility edge, labeled 
M.E., while the upper band edge remains at fj, — 2. 



plots the same as the dashed line but with a second lat- 
tice, V n , with V — 0.5 for chemical potentials near the 
upper mobility edge (labeled M.E.). < X > oa remains 
zero where we expect extended states but pins near A for 
fi > 2 — 2V. The mesoscopic version of the gapless insu- 
lator results in the early onset of pinning in the regime 
2 — 2V < fi < 2 and Eq. [3 Furthermore, the localized 
states with the additional lattice, V = 0.5, also display 
weak periodicity in < X > oc as a function of /x. These 
oscillations correspond to the chemical potential passing 
through peaks and valleys in the corrugated confinement 
potential. 

Fluctuations in the lattice depth can soften the other- 
wise sharp mobility edge. The quantity of interest, 2V/t, 
can fluctuate wildly with only moderate changes in Vl 
at extremely large lattice depths. To see this consider an 
approximate expression in terms of the hopping extracted 
from an analysis of the related Mathieu problem: V/(rj 
(V^Vy4)(V L /£^)- 3 / 4 exp(2 v /TV^)- A relative error 
in V and Vl, Ry and Ry L respectively, propagates to a 
relative error in 2V/t: [R\ + #^(3/4 - ^Vl/Er) 2 ] 1 / 2 . 
We have checked that this formula is quantitatively ac- 
curate for Vl ^ 5-Er by comparing with error derived 
numerically from the exact tunnelling. We find that for 
Rv = Rv L — 5% the relative error in 2V/t remains below 
20% for V L < 20E R . 

We note that additional time dependence in the model 
discussed here possesses other applications. We take Hcl 
as a good approximation to the nearly-periodic Harper 
model in the limit Eq. In the presence of a pulsed 
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secondary lattice: V oc J^ . 6 (T — jTo), where for inte- 
ger j the secondary lattice oscillates with period T , we 
simulate the kicked Harper model via Hql- The kicked 
Harper model exhibits chaotic behavior with the "classi- 
cal" to quantum crossover controlled by e. 

We have explicitly demonstrated the existence of a mo- 
bility edge (and the associated, unusual metal-insulator 
transition in a deterministic disorder-free environment) 
in suitably designed aperiodic cold atom optical lattice 
systems. The deterministic aperiodic background poten- 
tial in these optical lattices leads to exotic and nontriv- 
ial energy eigenstates dependent on the relationship be- 
tween irrational numbers and their rational approxima- 
tions. The ensuing quantum localization occurs in the 
absence of disorder and therefore distinguishes itself from 
Anderson localization which, in the solid state, masks the 
presence of mobility edges formed from quasiperiodic po- 
tentials in one dimension. 
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